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Abstract 

When the cancellation of the leading renormalon contributions is incorporated, the total 
energy of a bb system, E tot b g(r) = 2 m po ie,b + Vqcd(^ )> agrees well with the potentials used 
in phenomenological models for heavy quarkonia in the range 0.5 GeV" 1 < r < 3 GeV" 1 . 
We provide a connection between the conventional potential-model approaches to the 
quarkonium spectroscopy and the recent computation based on perturbative QCD. 



1 Introduction 



For over 20 years, most successful theoretical approaches for describing the charmonium and 
bottomonium systems (including the excited states) have been those based on various phe- 
nomenological potential models. These phenomenological-model approaches have elucidated 
nature of the heavy quarkonium systems, such as their leptonic widths and transitions among 
different levels, besides reproducing the energy levels. The phenomenological potentials deter- 
mined and used in these studies have more or less similar slopes in the range 0.5 GeV -1 r < 
5 GeV" 1 , which may be represented by a logarithmic potential oc logr + const. See e.g. Ref.|IJ 
for a recent analysis based on the potential models. An apparent deficit of these approaches is, 
however, a difficulty in relating phenomenological parameters to the fundamental parameters 
of QCD. 

The reason why people have been using phenomenological models is because the theory of 
non-relativistic boundstates based on perturbative QCD failed to reproduce the charmonium 
and bottomonium spectra. This is in contrast to the corresponding theory based on perturba- 
tive QED, which has been successful in describing the spectra of the QED boundstates. The 
main problem has been poor convergence of the perturbative expansions when the energy levels 
of the heavy quarkonia are computed in series expansions in the strong coupling constant. Since 
the coupling constant is quite large at relevant scales, approximating order one, it has been con- 
sidered as an indication of large non-perturbative effects inherent in these quarkonium systems. 
In fact the difference between a typical phenomenological potential and the Coulomb potential 
tends to be a linearly rising potential at distances r > 1 GeV -1 , suggesting confinement of 
quarks. Within perturbative QCD, the origin of the poor convergence has been understood in 
terms of the renormalon contributions ||. 

More recently, theoretical frameworks based on QCD have been developed for describing 
these quarkonium systems systematically. Within effective theories based on appropriate expan- 
sions in small parameters, various potentials are defined such that the leading-order potential 
plays a role close to that of the potentials introduced in the above phenomenological approaches. 
The order countings of terms in organizing the expansions depend crucially on the relative sizes 
of the dynamically generated scales (soft scale ~ mv , ultrasoft scale ~ mv 2 , where m and v are 
the quark mass and velocity, respectively) and the hadronization scale ~ Aqcd- This aspect 
contrasts with the fact that the expansion parameter in the non-relativistic boundstate theory 
based on perturbative QCD is simply 1/c (inverse of the speed of light). In the formalism 
developed in 0, f|, 0, @, 0|, or in the potential-Non-relativistic QCD (pNRQCD) formalism 
H, H formulated more systematically, potentials are defined in a way suited for practical com- 
putations by lattice simulations (or by using models). On the basis of these formalisms, lattice 
calculations have shown from first principles that the leading-order potential has a shape con- 
sistent with the phenomenological potentials in the relevant range, although the accuracy of 



the computations needs further improvements |1C, 11 



Very recently, a new computation of the charmonium and bottomonium spectra has been 
reported in the framework of non-relativistic boundstate theory based on perturbative QCD 
fl2| . It incorporated recent significant developments in the field: (1) the full computations of 
the quarkonium energy levels up to order 1/c 2 |13|, 111, ITR KB; (2) the cancellation of the leading 
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renormalons contained in the quark pole mass and the static QCD potential [T7|, IB]. As a result, 
convergence property of the series expansions of the energy levels improved drastically, which 
enabled stable perturbative predictions for the levels up to some of the n = 3 bottomonium 
states and the n = 1 charmonium states (n is the principal quantum number). Furthermore, the 
computed spectrum, when averaged over spins, reproduced the gross structure of the observed 
energy levels of the bottomonium states, within moderate theoretical uncertainties estimated 
from the next-to-leading renormalon contributions. It indicates that non-perturbative contri- 
butions to the bottomonium spectrum, in the scheme free from the leading renormalons, would 
absorb the next-to-leading renormalon uncertainties of the perturbative predictions and may 
be of the size comparable to them. 

It is then natural to ask whether there is a connection between the above phenomenologi- 
cal potential-model approaches (supplemented by the more systematic frameworks and lattice 
calculations) and the recent computation based on perturbative QCD. Once this connection is 
established, we may merge these approaches and further develop understandings of the char- 
monium and bottomonium systems. For instance, in the perturbative computation, the level 
splittings between the .S-wave and P-wave states as well as the fine splittings among the nPj 
states turn out to be smaller than the corresponding experimental values. Although the discrep- 
ancy is still smaller than the estimated theoretical uncertainties of the perturbative predictions, 
it should certainly be clarified whether they are explained by higher-order perturbative correc- 
tions, or, we need specific non-perturbative effects for describing them. On the other hand, 
the potential approaches have been successful also in explaining the S-P splittings and the fine 
splittings. Hence, we expect that a connection between these theoretical approaches would help 
to clarify origins of the differences of the present perturbative predictions and the experimental 
data. 

In this paper we focus on the perturbative static QCD potential up to (9(a|), since it dictates 
the major structures of the quarkonium spectra in the perturbative computation up to 0(1/ c 2 ) 
[ 12| . Taking into account the above key ingredient (2), we subtract the leading renormalon 
contribution from the QCD potential. Then we compare it with the phenomenologically deter- 
mined potentials. Our comparison also elucidates to which extent the perturbative computation 
of the QCD potential [up to C(a|), and after subtracting the leading renormalon] reproduces 
the results of the non-perturbative computations. (We will regard typical phenomenological po- 
tentials as representatives of the lattice results, taking into account consistency of the potentials 
determined in both approaches.) 

In Sec. 2 we review the theoretical uncertainties from the renormalon contributions within 
the context of the large-/9o approximation. In Sec. 3 we analyze the total energy of a quark- 
antiquark system up to 0(a%). Also the interquark force is analyzed in Sec. 4. We draw 
conclusions in Sec. 5. 



2 



2 Renormalons in the Large-/?o Approximation 



The static QCD potential, denned from an expectation value of the Wilson loop, represents the 
potential energy of a static quark-antiquark pair: 



VqcdM = — h m loj 



(0 Tr P exp 






0) 


(0 


Tr 1 


0) 



(1) 



= -C F av{1/r \ (2) 
r 

where T is a rectangular loop of spatial extent r and time extent T. The second line defines the 
F-scheme coupling constant, ay(l/r), where Cf = 4/3. In perturbative QCD, the V^-scheme 
coupling constant is calculable in a series expansion in the coupling constant asf] 

a v (l/r) = a 5 (/i) f) P n Qag(jir)) (^^j . (3) 

Throughout this paper, q s (/i) denotes the strong coupling constant in the MS scheme with m 
active flavors; \l is the renormalization scale; P n {L) denotes an n-th-degree polynomial of L. 
Although the exact QCD potential Vqcd(^) is independent of the scale //, at each order of the 
perturbative expansion /i-dependences remain. We keep fi as a free parameter in this section. 
From an analysis of higher-order terms, it has been known that the perturbative expansion 
of VqcdO") has an uncertainty of order Aqcd, which is referred to as the renormalon problem. 
We first review this property and estimate uncertainties of the perturbative prediction for the 
QCD potential. (See e.g. |2(| |2]J for introductory reviews.) 

The "large-/3o approximation" |22J is an empirically successful method for analyzing large- 
order behaviors of physical quantities in perturbative QCD and renormalon ambiguities inherent 
in them. Let us denote by Va (r) the QCD potential within this approximation and by V^\r) 
its 0{a n s +1 ) term: 

oo 

tfc(r) = £< } (0- (4) 

n=0 

From the Taylor expansion of the Borel transform of Vp Q (r) , we can easily compute (r) one 
by one from the lowest order. Also the asymptotic form for n ^> 1 is determined as 

V#»(r) ~ -C M ») x tg. x {^}" „!, (5) 

where (3$ = 11 — 2rii/3 is the coefficient of the QCD one- loop beta function. The above 
asymptotic behavior is independent of r. It means that, although each term of the potential is 



* From 0(ag) and beyond, the series includes infrared divergences; the divergences can be circumvented by 
a resummation of diagrams, which brings in log as in the series expansion, or, log(/x e ff r) term when the theory 
is matched to the pNRQCD effective theory [p~9|, p| . 
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Figure 1: The QCD potential in the large-/?o approximation truncated at 0(ctg +1 ) term. We set i± = 2.49 GeV, 

ni = 4 and as(n) = 0.273 [corresponding to a^(Mz) = 0.1181]. (a) Before subtraction of the leading 
renormalon. (b) After subtraction of the leading renormalon. 



a function of r, its dominant part for n ^> 1 is only a constant potential which mimics the role 
of the quark mass in the determination of the total energy of a quark- ant iquark system. As 
we raise n, first |Vg™ (r)| decreases due to powers of the small as', for very large n it increases 

due to the factorial n\. Around n = 27r/(/3 as(/i)), |Vg^(r)| becomes smallest. The size of 
the term scarcely changes within the range n G (n Q — ^/nQ,riQ + ^/no). We may consider the 
uncertainty of this asymptotic series as the sum of the terms within this range, since one may 
equally well truncate the series at order n — ^Juq or at order n + ^/hq in estimating the "true 
value" of the potential: 



SVfh(r) 



no- 



A 



/i exp 



2tt 



Poa s (n) 



(6) 



The /i-dependence vanishes in this sum, and this leads to the claimed uncertainty In Fig. |T](a) 
we show the QCD potential in the large- /3q approximation truncated at the (iV+l)-th term, 
J2n=o ^A^( r )' f° r N = 0,1,2,... and fi\ = 4. We see that the higher order corrections are 
indeed large and almost constant (independent of r). 

It was found |17], pl| that the leading renormalon contained in the QCD potential gets 



cancelled in the total energy of a static quark- ant iquark pair, 



Etot(r) = 2m po ie + V ( 



QCD l 



(7) 



if the pole mass m po i e is expressed in terms of the MS mass. Namely, when expressed in terms of 
the MS mass and in a series expansion in as(fi), the pole mass contains the leading renormalon 
|23| which is one half in size and opposite in sign of the leading renormalon of Vqcd(^)- Thus, 
the total energy E tot {r) is free from the leading renormalon uncertainties. E tot (r) possesses a 
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residual uncertainty originating from the next-to- leading renormalon [[|, 



6E tot (r) ~ A x (Ar) 2 , (8) 

which is smaller than the leading renormalon uncertainty in the range r < A -1 . Shown in 
Fig. |l](b) is the QCD potential in the large-/3o approximation [truncated at the (AT+l)-th term] 
after the leading renormalon is subtracted at each order of «s(A i ) : 



n=0 



Vp ( r ) + °f ^asiv) ■ • { 2?r ) n\ 



(9) 



One sees that the series expansion of the potential has become much more convergent as 
compared to Fig. |l](a). For a particular choice of the scale /i = 2.49 GeV, the term on the right- 
hand-side of Eq. (0) becomes smallest at around n = 7 in the range 1 GeV" 1 < r < 5 GeV" 1 . 
Hence, the error bars corresponding to the next-to-leading renormalon uncertainty ±|A • (Ar) 2 
(taking A = 300 MeV) are attached to the potential for = 7 in the same figure. We may 
consider that the line for N = 7 together with the error bars indicate a typical accuracy of the 
perturbative prediction for the QCD potential, when the leading renormalon is cancelled. We 
see that the potential is bent upwards at long distances as compared to the leading Coulomb 
potential (N = 0). If we choose a smaller scale for /i, the term becomes smallest at a smaller 
n. In this case, convergence properties become better at larger r, where we obtain a value of 
V j3 (r) consistent with N = 7 of Fig. 0(b) with less terms (smaller N) . Similarly to the leading 
renormalon case, the uncertainty is ^-independent, nonetheless. 



3 The Total Energy of a qq System 



Now we examine the total energy of a quark-antiquark pair, defined in Eq. ([?]), exactly up to 
0{ot z s ). This quantity is free from the leading renormalon uncertainty; in fact the cancellation 
of the leading renormalons occurs at a deeper level than what can be seen in the large-/3o 
approximation ||18|| . We also note that the cancellation at each order of perturbative expansion 



is realized only when we use the same coupling constant in expanding m po i e and Vqcd(^)-Q 
The QCD potential of the theory with n/ massless flavors onlyQ is given, up to (9(a|), by 



V ( 



QCD l 



a 




+ 2(/3i + 2/5 ai)£ + a 2 



(10) 



t This can be seen, for example, from the fact that the order no = 2-k / ((3oas(n)) at which Eq. (|^) becomes 
smallest is dependent on the value of as(p) used for the expansion. 

* The QCD potential of the theory which contains rih heavy flavors (with mass m) and ni massless flavors 
coincides with the potential in Eq. jl0| ) up to 0{a 3 s ) if we count 1/r = O(asm) and if we properly match the 
coupling to that of the theory with ni massless flavors only. 
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where PI 



Z = \og{iir) + 1 E, (11) 

A> = 11 - ^n,, A = 102 - yn h (12) 

31 10 4343 2 9vr 4 /1229 52 Cs\ 100 2 

«i = y-y^, a 2 = — + 36. +66(3-^-^ ^_ + _j+_ ni . 

(13) 

The relation between the pole mass and the MS mass has been computed up to three loops in 
a full theory, which contains nn heavy flavors and n\ massless flavors [p5fl . (The same relation 
was obtained numerically in in a certain approximation.) Rewriting the relation in terms 
of the coupling of the theory with n\ massless flavors only, we 




m pole = m \ 1 + \ + ( I di + I I d 2 V , f 1 I ) 



3 7T 

where m = ^Ms( m Ms) denotes the renormalization-group-invariant MS mass, and 

3049 2tt 2 vr 2 log2 s> . ,. „ . . , „ . 

dx = - W + -g- + —f--¥ + n l [-—--)+n h ( — + -), (15) 

1145453 25379 tt 2 

d 2 = 1 h 

10368 2592 

19 log 4 2 76 a 4 




54 9 432 216 



/ 81227 965tt 2 ll7r 2 log2 707Cs 61 tt 4 2vr 2 log 2 2 log 4 2 8a 4 
+Hl \ 7776 ~ ~648 81 216~ + 1944 + 81 + 81 + ~27~ 

/ 157007 13627 vr 2 640 tt 2 log 2 751 Cs 41 vr 4 vr 2 log 2 2 log 4 2 

+Uh 1 1 1 1 

I 7776 1944 81 216 972 81 81 



8 a 4 
+ 27 " 
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5C 
4 


2353 
23328 + 


13 TT 2 

324 


+ 7 C 3 

54 


9481 


4tt 2 


11 Cs 


23328 


405 


54 



+ n l I „oo^o + -r^rr + I + n i n h 



( 5917 13 tt 2 2Cs 



I" 



v 11664 324 27 

w im. _ u! _ ii^S ' (16) 

fc \23328 405 54 /' v ; 

with a 4 = Li 4 (|). Furthermore, we rewrite as(m) in terms of a-sil 1 ) using the renormalization- 
group evolution of the coupling constant. Thus, we examine the series expansion of E tot (r; m, cts(/i)) 
in cxs{fj) up to 0(oP s ). Qualitatively the series shows a convergence property very similar to 
VpJr) foriV = 0,l,2; cf. Fig. |(b). 

§ When rih — 1, this relation coincides with Eq.(14) of p5fl , which is given numerically (indirectly through 
/?o and /3i). Note that, in the other formulas of |p5[ , the coupling of the full theory is used. 
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The obtained total energy depends on the scale \i due to truncation of the series at a finite 
order. One finds that, when r is small, the series converges better and the value of E tot (r) is 
less /i-dependent if we choose a large scale for /i, whereas when r is larger, the series converges 
better and the value of E tot (r) is less /i-dependent if we choose a smaller scale for /i. Taking 
into account this property, we will fix the scale [i in two different ways below: 

1. We fix the scale fi = /^i(r) by demanding stability against variation of the scale: 

= 0. (17) 



H-^-E tot (r;m,a s (lJ,)) 
djj, 



2. We fix the scale \i 
[£>(a|) term] of E, 



/i2(r) on the minimum of the absolute value of the last known term 



tot l 



d_ 



fJ.=V2 (r) 



In this analysis we examine the total energy of a bb system. We set mj, = m^ /[S (m^ IS ) = 
4.203 GeV, which is taken from (Its error is estimated to be about ±30 MeV.) For 

simplicity we analyze E tot {r) in two hypothetical cases: (i) when m c = (jii = 4 and rih = 1), 
and (ii) in the limit m c — > mj (ni = 3 and = 2). The real world lies somewhere in between 
the two cases: the charm quark decouples in the excited states of bottomonium but not in 
the ground state ||27|| . A more precise analysis requires inclusion of nonzero m c effects into 
Etot{ r ), which will be reported elsewhere [p8 |. The input value of the strong coupling constant 



is as\Mz) = 0.1181 f29|. We evolve the coupling and match it to the couplings of the theory 
with ni = 4 and 3 successively by solving the renormalization-group equation numerically with 
the 3- loop beta function and by using the 3-loop matching condition (3- loop running). 

Fig. |^(a) shows E tot {r) (measured from 2Wib) for the two cases (i) and (ii). In each case 
-E'tot(^) are plotted with the two different scale- fixing prescriptions; the total energy hardly 
changes whether we choose /i = /ii(r) or fi = ^{f)- In case (i), the minimal sensitivity scale 
//i(r) exists only in the range r < 3 GeV -1 ; for the choice = /i2(7"), the minimum value of 
\Etot( r )\ is zero i n the range r < 3 GeV -1 , whereas \E^l(r)\ > in the range r > 3 GeV -1 . 
These features indicate an instability of the perturbative prediction for E tot {r) at r > 3 GeV -1 . 
The scales Hi(r) and ^{r) are shown as functions of r in Fig. H(b). For comparison, we also show 
/i = exp(— 7_e)/r, which has been considered as a natural scale of the QCD potential, Vqcd(^), 
conventionally. One sees that /ii(r) and ^{r) are considerably larger than exp(— 7^)/r. The 
scales chosen in case (ii) are similar. In Table [l| we show each term of the series expansion of 
E tot (r). The series shows healthy convergent behavior at r < 3 GeV -1 . 

At this stage, let us discuss why the scales /ii(r) and /i2( r ) are considerably larger than 
exp(— 7e)/r. For this purpose we use an approximate expression for the pole mass, which 

™ We take the matching scales as rf%b and m c (= mi,), respectively. 



7 



4 




Figure 2: (a) The total energy of a bb system measured from 2mf, in two hypothetical cases. In each case, the 
scale is fixed by \x — [i\ir) (dotted lines) or \i — (solid lines if E to ' t (r) = 0; dashed lines if |£^ t(r)| > 0). 

(b) The scales chosen by the scale- fixing prescriptions (|l7]) and ( |ilf ) in case (i). The notations are same as in 
(a). A conventional scale choice /i = exp(— 7s)/r is also shown (dotdashcd line). 





case (i) 


n = m(r) 


H = /i 2 (r) 




E£l(r) 


E£l(r) 






E£l(r) 


r = 1 GeV" 1 


797 


69 


-17 


750 


98 





r = 2 GeV" 1 


1255 


-14 


-17 


1173 


48 





r = 3 GeV" 1 


1709 


-290 


13 


1606 


-185 
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case (ii) 


[i = n x [r) 


H = n 2 {r) 








E^(r) 






r = 1 GeV" 1 


962 


70 


-32 


879 


117 





r = 2 GeV" 1 


1659 


-116 


-31 


1502 


4 





r = 3 GeV~ x 








1994 


-197 


70 



Table 1: Series expansion of the total energy in as{(J-) with the two scale choices Eqs. (|l7| ) and ([FJ|). e[^1{t) 
denotes the C(a§) term of E tot (r). All numbers are in MeV unit. The minimal sensitivity scale Hi(r) exists 
only at r < 2.8 GeV -1 in case (ii). 
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follows from the fact that the dominant contribution to the pole-MS mass relation can be read 
from the infrared region, loop momenta g C m, of the QCD static potential Hl8|: 



2m 



pole 



/d^q 2CJp f m 

|Vqcd(<?)| = 2m + — — j dqa v (q). 



2C F 



(27T) 



(19) 



|cj|<m 



Here, Vqcd(?) — —Cp47iav(q)/q 2 is the QCD static potential in momentum space. Then the 
total energy can be written approximately as 



Etot(r) 



f d^q r 
2m + / ,„ NO |VQ C D(g)| 0(m - q) - exp(iq ■ r) 



(2tt)* 



9- m 2Cf 
2m H 



7T Jo 



dqay(q) 



9(m — q) — 



sin(gr) 



qr 



(20) 
(21) 



In the integrands, the factors in the brackets [• ■ •] are appreciable only in the range 1/r < q < m. 
So, roughly speaking, E tot (r) is determined from an average (ay) of the ^/-scheme coupling 
cty(q) over the range 1/r < q < m. When evaluating this quantity in fixed-order perturbation 
theory, a scale fi{r) which represents this average coupling, i.e. oy(//(r)) ~ (ay), would be a 
most natural scale. Such a scale should lie between 1/r and m. This argument is in contrast 
with the conventional principle for the scale choice for the QCD potential Vqcd(^)- Apart 
from Aqcd, the QCD potential contains only one scale 1/r, so that the choice of scale has 
been almost automatic, [i ~ 1/r. The potential alone, however, has a large uncertainty due to 
the leading renormalon. It stems from the contribution of &v{q) at q ~ Aqcd- On the other 
hand, the total energy is free from the leading renormalons by cutting out large contributions 
from Aqcd ~ g < 1/r as seen in Eq. (pi)) . Consequently the relevant scale is shifted to higher 
momentum region in comparison to that of Vqcd(?")- 

It would also be instructive to compare the above scale choices with the Brodsky-Lepage- 
Mackenzie (BLM) scale-fixing prescription |30| applied to Vqcd(^) an d -Etot( r ), respectively. In 
this prescription (at the lowest order), the part of higher-order corrections to VqcdM or to 
E tot {r) given by the large- (3 approximation is absorbed into the scale choice. For the QCD 
potential, at the lowest order the BLM scale is fixed as \i = exp(— |— 7e)/r ~ 0.43 exp(— 7s)/r. 
For the total energy, the BLM scale at the lowest order is given by /i = f(mr)/r, where 



f(x) = exp 



(logx 



7[- 



+ IE 



53 



V 



x 



192' x 



5 
6 



IE 



[22) 



Due to the singularity of f(x) at x = ir/2, the BLM scale turns out to be unstable around 
r = 7r/(2m). This is because the coefficient of /3 log// in E tot (r;rn,as{^)) becomes small by a 
cancellation between VqcdM an d 2m po i e . In this region of r, the BLM prescription for E tot (r) 
would be unreliable. For r > 7r/(2m), the function f(mr) increases monotonically. Setting 
fn = 4.203 GeV, we find that at r > 0.6 GeV -1 the scale fi = f(mr)/r exceeds the BLM 
scale of the QCD potential; at r > 1 GeV -1 , [i = f(mr)/r becomes almost independent of 
r, fi ~ 0.5 GeV, converging towards the BLM scale of the pole mass. These features in the 
region r > 0.6 GeV -1 are consistent with the results of the analysis given in Sec. |2|. Since the 
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Power-law Pot. 



1.5 



> 

o 1 



0.5 



Logarithmic Pot. 
£ tot (r):case (ii) 




[GeV" 1 ] 



¥(1S) 



Figure 3: A comparison of the total energy of a bb system in the two hypothetical cases (dotted lines) and 
typical phenomcnological potentials (solid lines). For a reference, we show typical sizes of the bottomonium and 
charmonium S states as determined from the r.m.s. interquark distances with respect to the Cornell potential: 



Cornell ' 



higher-order corrections are large for Vp (r), the BLM scale of Vqcd(^) tends to be small. On 
the other hand, since the higher-order corrections are smaller for V/3 (r), the BLM scale of 
Etot{r) is larger. Intuitively the BLM scale of a quantity sensitive to renormalons is attracted 
towards the Aqcd scale, whereas that of a renormalon-free quantity is determined by a short- 
distance scale. Thus, the qualitative features of the BLM scales agree with those of the scales 
shown in Fig. 0(b) in the range 1 GeV -1 ^ r < 5 GeV -1 , although the level of agreement is 
not very accurate. At shorter distances r « 7r/(2m), validity of the BLM prescription for the 
total energy seems doubtful on account of the large cancellation. 

We return to the discussion of E tot (r). In Fig. [3] we compare the total energies in case (i) 
and (ii) with typical phenomenological potentials used in phenomenological approaches. We 
take: 

• A Coulomb-plus-linear potential (Cornell potential) |32| : 

V(r) = -- + - 2 (23) 
r er 

with k = 0.52 and a = 2.34 GeV -1 . 

• A power-law potential |j33"| : 

V{r) = -8.064 GeV + (6.898 GeV)(r x 1 GeV) 01 . (24) 
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A logarithmic potential EM : 



V(r) = -0.6635 GeV + (0.733 GeV) log(r x 1 GeV). (25) 

We may consider the differences of these potentials in the range 0.5 GeV" 1 < r < 5 GeV" 1 
as uncertainties of the phenomenologically determined potentials. In order to make a clear 
comparison, arbitrary constants have been added to all the potentials and E tot (r) such that 
their values coincide at r = 1 GeV" 1 . As stated, we expect the perturbative prediction for 
a realistic E tot {r) to lie between those for the cases (i) and (ii). It appears to be in good 
agreement with the phenomenological potentials in the above range. The level of agreement is 
consistent with the uncertainties expected from the next-to-leading renormalon contributions 
(indicated by the error bars). 

4 The Interquark Force 

Instead of the total energy, we may also consider the interquark force defined by 

F{r) = —E^ir) = -^V QC n(r) (26) 

= -C F aF ^ r \ (27) 

The last line defines the "F-scheme" coupling constant a F (fi). The interquark force is also free 
from the leading renormalon. In fact it has been noted ]35| that the perturbative expansion of 
F(r) is more convergent than that of the potential Vqcd(^)- Since in the zero quark mass limit 
F(r) is dependent only on r, we may determine its r-dependence using the renormalization- 
group equation: 

l^-^apin) = f3p{otp). (28) 

It is instructive to compare the beta functions for the couplings defined in the three different 
schemes. For n\ = 4, we find 

Pv(a v ) = - 0.6631 a 2 v - 0.3251 afy - 1.7527 a£ + 0(a 5 v ) (V-scheme) 
(3 F {a F ) = - 0.6631 a£ - 0.3251 a% - 0.5861 a% + 0(a%) (F-scheme) (29) 

p m (a s ) = - 0.6631 a| - 0.3251 a| - 0.2048 a| + 0{a%) (MS-scheme) 

The first two coefficients of the beta functions are scheme-independent (when we neglect the 
quark masses). The third coefficient of the ^/-scheme beta function is quite large, reflecting poor 
convergence of VqcdM due to the leading renormalons. The third coefficient of the F-scheme 
beta function is smaller by factor 3 due to cancellation of the leading renormalon. The third 
coefficient of the MS-scheme beta function is even smaller by factor 3. This may be due to the 
fact that the F-scheme coupling still contains the next-to-leading renormalon contributions. 
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From this comparison, we may conclude that it is better to analyze F(r) rather than VqcdO") 
physical quantity, in perturbative analyses.Fl 

The observed bottomonium spectrum is qualitatively very different from the Coulomb spec- 
trum. The largest difference is that, the level spacings between consecutive bottomonium nS 
states are almost constant, whereas in the Coulomb spectrum the level spacings decrease as 
1/n 2 . When we consider effects of the QCD radiative corrections on the lowest-order Coulomb 
potential, one may interpret that in the QCD potential, —Cpay(l/r)/r, the ^-scheme coupling 
increases at long distances, so that the potential will be bent downwards. This is obviously a 
bad interpretation, because in such a case, the level spacings among the excited states become 
even smaller than those of the Coulomb spectrum. We should rather consider the interquark 
force. A better interpretation is that in F(r) = — Cpap(l/r)/r 2 , the F-scheme coupling in- 
creases at long distances, and |F(r)| grows correspondingly. This means that the slope of the 
potential becomes steeper at long distances. (Its effect resembles an addition of a linearly rising 
potential to the Coulomb potential.) Accordingly the level spacings among the excited states in- 
crease. Thus, the effects of the radiative corrections on the level spacings are even qualitatively 
reversed, whether we consider VqcdO") ot F(r) as the physically relevant quantity.fj] 

One may verify these features in Fig. |], in which the Coulomb potential, the ^-scheme 
potentials and the F-scheme potentials are displayed. The ^-scheme potentials are calculated 
by solving the renormalization-group equation for ay numerically, using (3y in Eq. ( |29| ) up 
to order a v (1-loop), order a v (2-loop) and order a v (3-loop). The F-scheme potentials are 
calculated by first solving the renormalization-group equation for ap numerically via j3p in 
Eq. (|29|) and then by integrating —F(r) over r numerically; arbitrary constants are added 
such that the F-scheme potentials coincide the Coulomb potential at r = 0.4 GeV -1 . The 
initial values for ay and ap are given at r = exp(— •yp)/ffib by matching to the fixed-order 
results. As can be seen, the ^-scheme potentials become singular at fairly short distances, 
r ~ 2 GeV -1 (1-loop), 0.9 GeV -1 (2-loop), and 0.4 GeV" 1 (3-loop), respectively. As expected, 
the F-scheme potentials have wider ranges of validity^]: they become singular at r ~ 6.9 GeV -1 
(1-loop), 2.8 GeV -1 (2-loop), and 1.7 GeV -1 (3-loop), respectively. The situation is puzzling, 
however, in that the predictable range reduces as we include more terms of /3p(ap). The 
2-loop and 3-loop F-scheme potentials are consistent with the phenomenological potentials 
within the uncertainty expected from the next-to-leading renormalon contributions, in the 
range 0.5 GeV" 1 < r < 2.8 GeV" 1 and 0.5 GeV" 1 < r < 1.7 GeV" 1 , respectively. On the 
other hand, the 1-loop F-scheme potential does not satisfy this criterion. 

If we take a larger input value for a ( g\Mz), the slopes of the F-scheme potentials get 
steeper, since ap increases. Also, it explains why F to t(r) for case (ii) is steeper than that for 

* This is valid up to the constant term of the potential, which is important in relating the boundstate masses 
to the heavy quark masses. An alternative way may be to study renormalization-group evolution of VqcdM 
after subtracting the leading renormalon from it by hand [similar to V(3 (r) of Eq. (0)]. 

t It is a matter of interpretaion. One may understand the radiative corrections in the context of ^-scheme 
and require for large non-perturbative corrections to remedy the discrepancy from the phenomcnologically 
determined potentials or the results of non-perturbative (lattice) calculations (see e.g. ]3(| |37j). Alternatively one 
may understand the radiative corrections in the context of F-scheme and call for much smaller non-perturbative 
corrections. 

t The large discrepancy between the potentials obtained from (3v{ay) and (3f(&f) was noted first in [B8|. 
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F- scheme Pot. 




V- scheme Pot. 



r [GeV 1 ] 



Figure 4: A comparison of the QCD potentials calculated in ^-scheme and in F-scheme as well as the Coulomb 
potential. The Coulomb potential is given by —Cpa/r with a — 0.279. The l^-scheme and ^-scheme potentials 
correspond to a^ 5) (Af z ) = 0.1181. 



case (i) in Fig. @(a): ap for ni = 3 is larger than that for ni = 4 at r > 1/mj,. 



5 Conclusions 



When we incorporate the cancellation of the leading renormalon contributions, the perturba- 
tive expansion of the total energy E tot {r) of a bb system, up to and supplemented by 

the scale-fixing prescription or (|18|), converges well at r < 3 GeV^ 1 . Moreover, it agrees 
with the phenomenologically determined potentials in the range 0.5 GeV -1 < r < 3 GeV -1 
within the uncertainty expected from the next-to-leading renormalon contributions. Even at 
r > 3 GeV -1 , the scale-fixing prescription (|l8| ) gives a reasonable prediction for E tot (r); it ap- 
pears that the perturbative prediction does not break down suddenly but rather the uncertainty 
grows gradually as r increases. The agreement is unlikely to be accidental, since as soon as we 
take the input a^\M z ) outside of the present world average values 0.1181 ± 0.0020 PD[ , the 
agreement is lost quickly. 

A non-relativistic Hamiltonian 



H = 2m po i e + p 2 /m po i c + V QC d(V) = p 2 /m poXe + E t , 



(30) 



constitutes a part of the full Hamiltonian (up to the order 1/c 2 ) analyzed in [12|. It detetermines 



the bulk of the quarkonium level structure computed therein. At the same time, the above 
Hamiltonian is exactly the ones analyzed in the conventional phenomenological potential-model 
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approaches (at the leading order) if E tot {r) is identified with the phenomenological potentials. 
Also, it is the leading-order Hamiltonian of the more systematic frameworks discussed in Sec. |I[ 
Thus, we find that the agreement of E^ ot {r) and the phenomenologically determined potentials 
is the reason why the gross structure of the bottomonium spectrum is reproduced well by the 
computation based on perturbative QCD. Our observation confirms the conclusion of [|l^], that 
once the leading renormalon contributions are cancelled, there remain no large non-perturbative 
effects, which essentially deteriorate perturbative treatment of some of the bottomonium and 
charmoninum states, but only moderate contributions comparable in size with the next-to- 
leading renormalons. 

Similarly, if we analyze the interquark force F(r) instead of Vqcd(^), the range of pertur- 
bative predictability becomes significantly wider, as known from the previous study p5| . We 
confirm this observation using a renormalization-group analysis. We find that the 2-loop and 3- 
loop renormalization-group-improved potentials, obtained by integrating —F(r), are consistent 
with the phenomenological potentials up to r ~ 2.8 GeV™ 1 and r ~ 1.7 GeV -1 , respectively. 

We expect that the connection elucidated in this work will be useful for developing deeper 
theoretical understandings of the bottomonium and charmonium systems. For more detailed 
comparisons, in general it would be more secure to compute the quarkonium spectra directly 
rather than E tot (r) or F(r). Indeed, the series expansions of the quarkonium energy levels turn 
out to be more convergent when we include the full corrections (p 4 -term, Darwin potential, 
spin-dependent potentials, etc.) to the 0{l/c 2 ) Hamiltonian, as compared to the expansions 
of the energy levels of the simplified Hamiltonian fl3!]| ) (even after the leading renormalons are 
cancelled) . 
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